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The  geometrical  optics  approximation  of  solutions  of  the  two  dimensional 
wave  equation  for  an  isotropic  inhomogeneous  medium  is  considered.  An  ordinary 
differential  equation  for  a  quantity  which  is  inversely  proportional  to  the 
geometrical  optics  field  intensity  (the  square  of  the  field  amplitude)  is 
derived.  This  new  equation  along  with  the  standard  ray  and  phase  equations 
form  a  system  from  which  a  complete  wave  solution  in  the  optical  limit  can  be 
calculated  numerically,  e.g.,  by  means  of  an  analog  computer  acting  as  a  differ- 
ential analyzer. 

The  examples  of  a  homogeneous  medium  and  a  medium  whose  index  of  refraction 
is  a  function  of  one  variable  are  discussed,  and  the  results  of  the  preceding 
analysis  are  verified  in  these  two  special  cases. 
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1,  Introduction 

The  method  of  ray  tracing  or  geometrical  optics  is  a  useful  tool  for  the 
theoretical  investigation  of  electromagnetic  wave  propagation  in  a  general 
inhomogeneous  medium.   If  the  propagation  takes  place  in  the  earth's  atmosphere, 
for  which  the  index  of  refraction  varies  quite  slowly,  the  method  ia  effective 
even  at  radio  frequencies.  By  the  use  of  analog  computers  in  connection  with 
ray  tracing  an  analysis  of  a  number  of  interesting  phenomena  observed  experi- 
mentally in  tropospheric  propagation  can  be  carried  out  on  models  for  which  a 
realistic  variation  of  the  index  of  refraction  may  be  assximed.*-  -• 

The  direct  information  one  can  obtain  by  the  use  of  the  ray  tracing  tech- 
nique is  in  the  form  of  the  ray  paths  and  the  approximate  phase  of  a  propagating 
wave.  If  these  quantities  can  be  obtained  analytically,  in  the  form  of  an 

explicit  solution  of  the  ray  equations,  one  can  also  obtain  the  approximate 

[21 
field  amplitude  analytically,'--'   If,  however,  nxunerical  integration  methods 

are  used  to  obtain  the  ray  paths  and  phase  the  field  amplitudes  are  less  easy 
to  find.  A  popular  way  of  estimating  the  relative  variation  of  field  amplitude 
in  space  is  to  observe  from  a  graphical  representation  the  spatial  variation  of 
the  density  of  the  rays.  In  addition  to  being  rather  inelegant  this  method  seems 
likely  to  be  less  accurate  than  is  necessary  since  there  does  exist  an  analytical 
connection  between  the  approximate  field  amplitude,  the  ray  paths  and  the  approx- 
imate phase. 

In  addition,  in  a  region  where  the  rays  become  extraordinarily  dense,  such 
as  the  neighborhood  of  a  caustic,  or  where  the  rays  disappear  entirely  in  skip 
zones  or  shadow  regions  the  graphical  procedure  for  obtaining  the  approximate 
field  amplitude  breaks  down. 


In  view  of  these  difficulties  it  seems  desirable,  therefore,  to  investigate 
the  possibility  of  extending  the  gecsnetrical  optics  method  for  treating  propa- 
gation in  an  inhomogeneous  medium,  especially  in  connection  with  numerical 
techniques  and  the  use  of  an  analog  computer.  This  article  will  deal  with  the 
problem  of  obtaining  the  approximate  field  amplitude  in  the  case  where  the 
field  of  rays  is  regular.  An  ordinary  differential  equation  which  is  coupled 
to  the  ray  equations  and  whose  solution  is  inversely  proportional  to  the  square 
of  the  field  ^plitude  in  the  geometrical  optics  limit  will  be  derived.  This 
differential  equation  can  be  integrated  numerically  by  means  of  the  same  proce- 
dures which  are  applied  to  the  ray  equations.  It  may  also  prove  to  be  of  some 
theoretical  value. 

A  future  article  will  deal  with  the  problem  of  extending  the  geometrical 
optics  theory  to  caustics  and  shadow  regions. 

This  article  viill  be  concerned  only  with  the  scalar  wave  equation  in  two 
dimensions.  In  the  optical  limit  polarization  effects  are  usually  negligible, 
and  for  practical  applications  the  two  dimensional  problem  is  often  sufficient. 
Considerable  difficulties  arise  if  an  attempt  is  made  to  apply  the  analysis 
given  here  to  the  three  dimensional  problem. 

For  the  sake  of  completeness  Section  2  of  this  article  will  deal  with  pro- 
perties of  the  geanetrical  optics  solution  of  the  wave  equation  that  are  more 
or  less  well  known.  In  Section  3  the  differential  equation  whose  solution  is 
inversely  proportional  to  the  square  of  the  geometrical  optics  field  amplitude 
is  derived.  In  Section  h   two  examples  are  given  to  illustrate  the  validity  of 
the  differential  equation. 


2.  The  Geometrical  Optics  Solution  of  the  Wave  Equation, 

The  wave  equation  in  two  dimensions  for  an  inhomogeneous  medium  whose 
optical  properties  are  given  by  the  index  of  refraction  n(x,y)  has  the  form 

(1)  u   +  u   +  k^n^u  -  0  , 

30C    yy  ' 


of  light  in  a  vacuum.   It  will  be  assumed  that  n(x, y)  is  positive  everywhere. 

The  standard  formal  procedure  for  obtaining  the  optical  approximation  for 
u  in  (1)  is  to  assume  a  form 

(2)  u(x,y)  -  w(x,y)e^^^  ^^'^^ 

and  to  substitute    (2)  into  (1).     The  result  Is 


(3)  k^in^  -    (Vt)^(w  +  ik^2Vt-  Vw  +  (Ai)wi>   +  Aw  -  0 


The  function  f  ,   called  the  elconal,  is  now  chosen  so  that  it  satisfies  the 
eiconal  equation 

(U)  (Vt)2  -  n^  . 

The  curves  ^  "  constant  are  the  geometrical  optics  wavefronts,  the  orthogonal 
trajectories  of  which  are  the  rays  of  the  optical  system.  Next,  it  is  assumed 
that  there  Is  a  function  A(x,y)  such  that  for  x  and  y  fixed 

is)  11m  w  •  A   and   lim  vw  =  VA  . 

k  -»a3  k  -XX) 

It  is  also  assumed  that 


(6)  lim     Aw   - 

k  ->oo    TT  ■  °  • 


-  u  - 

After  the  application  of  (h),    (5)  and  (6)  to  (3)  there  results 

(7)  2Vt.  VA  +  (At)  A  =  0  , 

which  is  the  transport  equation  for  the  field  amplitude  A  in  the  optical  limit, 

The  first  order  partial  differential  equation  (Li)  for  the  eiconal  is 
equivalent  to  a  system  of  ordinary  differential  equations  for  the  rays  of  the 
optical  system  according  to  the  theory  of  Hamilton  and  Jacobi.     The  ray 
equations  are 

X  =  p  ,  p  -  ^(n  )^ 

(8)  y  "  q  ,  q  -  kn^y 


2^2    2 
p  +  q  »  n 


where  the  dots  indicate  differentiation  with  respect  to  a  parameter  t  in  terms 
of  which  the  rays  are  described: 

X  -  x(t  ) 

y  =  y(T)  . 

The  last  equation  of  (8)  is  an  integral  of  the  first  four  equations,  so  that 

it  holds  for  all  values  of  t  if  the  initial  values  of  the  system  at  some  value 

T  of  T  are  selected  in  such  a  way  that  it  holds  initially: 
o 


p2(T,)  *q2(T^)  =n2[x(T^),  y(.^)]  . 


For  t^  =  (n^)^  X  +  (n^)  y  »  2(pp  +  qq)  -  ^  (p^  +  q^) 


-  5  - 

The  eiconal  ^  {x,y)   is  thus  given  by 


(9) 


T 

t(x,y)  -  J  n^[x(T),y(T)]dT  +  t  [x(T^),y(T^)]  , 


.o 

a  fact  which  can  be  demonstrated  by  showing  that 

(10)  P  ■  V  ^  ""^y  • 

It  is  convenient  for  the  purposes  of  this  article  to  use  the  eiconal  \^, 
itself,  rather  than  t  as  the  ray  parameter.  To  do  this  let  us  consider  the  sys- 
tem generated  by  the  prescription  of  initial  values  on  some  initial  curve 

(11)  +(x,y)  *  t  =  constant  , 

which  is  an  initial  value  of  the  eiconal,   and  describes  the  initial  position  of 
a  wavefront.     We  can  suppose  that  the  curve  described  by   (11)   can  be  expressed 
in  terras  of  a  parameter  C: 

x(tQ,o-)  -  xjc^)  P(V<5-)  =  p^(oO 

(12)  with  and  p^(o')+q^(or)=n^(x^,y^)   . 

y(t^,cr)  ■=  y^(cr)  q{f^,<s-)  -  q^(cr) 

The  ray  solution  of  (8)  will  be  a  pair  of  functions  x(i,cr),   y(t,o')  having 
the  property  that 

(13)  x'x^i,    y  -  y^  ^,  p  -  p^if,  q  -  q^t 

since  t  and  CT  will  be  independent  variables.  From  (9)  and  (13)  the  equations 
(8)  can  be  written 


OR 


WAVEFRONT  xj/  =  i/^^ 
X  =  XJo") 


INITIAL    WAVEFRONT    SURFACE    i// =  i//^ 


6  - 


2n 


(n^l 


'♦'    n     ^    2n 


2    2    2 
p  +  q  =  n  . 

It  can  be  demonstrated  easily  that  vdth  the  initial  conditions  (12)  and 
if  the  first  four  equations  of  (1)4)  are  satisfied,  the  last  equation  of  (Hi) 
holds  automatically.  For  from  the  first  four  equations  of  (U)  it  follows 
that 

(15,         ^^.  5^.^.  ^,,  -  .p^  ,  a,^  .  ^  ,p=  .  ,h  . 

After  an  integration  of   (15)  with   respect  to  iff    there  results 

(16)  p2  +  q2  -  n^   +  f  (c^)   , 

where  f (C)  is  a  function  of  C  alone.  But  if  t  is  set  equal  to  t   in  (16),  then 
the  last  equation  of  (Hi)  implies  that  f{d)   -  0,  and  the  demonstration  is  com- 
pleted. 

It  can  be  demonstrated  similarly  with  the  aid  of  (Ih)  that  the  vector 
(p,q)  is  orthogonal  to  the  vector,  (x,y,)  everywhere  if  this  is  true  initially. 
Consider 


g|  (px^  *  qy^)   -  Pf,x^  -H  q^y^  ^  p^^  -.  qy,^^^ 
(17)  .lli..^.l!!iz.    .,4/4Vq4J4 


r  ^o  "^  p  3^  -  r  ^  a?  T 


2n'^       "         2n"       °  °"  V"    ; 


-     1     3(n^)   ^     1     a(n^)         1  a(n^)  .   _ 

rr-do^*  -2-5^-  J~^5T'  °  * 

2n  2n  n 


where  again  (lU)  has  been  used.   It  follows  from  (1?)  that 

(3B)  px^  +  qy^  »  g(cr)  , 

where  g(a)   is  a  function  of  o  alone.  If  the  initial  conditions  in  (12)  are 
chosen  so  that 

(19)  Po^^^^^o^''^^  *  %(^)yo^<=^) "  ° » 

where  the  primes  indicate  differentiation  with  respect  to  o",  then  by  setting 

f  =   ^    ±n   (19)  it  follows  that  g(cr)  -  0:  this  completes  the  demonstration.  This 
o 

orthogonality  condition  as  well  as  the  last  equation  of  (12 )  can  be  satisfied 


initially  by  the  choice 


(20)  [  ]-  I 


Now  the  vector  (t  ,t  )  is  clearly  orthogonal  to  the  vector  (x^y^  every- 
where, for 

(21)  °-|?-^x  V^^V 

Hence  the  vector  (p,q)  is  proportional  to  the  vector  (t  ,t  ).  Moreover,  from 
(lU)  there  results 

pijr   +  q^ 

(22)  —^ ^  'K^f'^y^    -  ^  • 

From  (22)  and  the  last  equation  of  (lU)  it  follows  that  the  proportionality 

factor  between  (p,q)  and  (t  .\|r  )  must  be  lh.9   constant  1.  Thus  the  relations  (10) 
A  y 


8  - 


hold  everywhere,  and  because  of  the  last  equation  in  (8)  it  follows  that  the 
eiconal  equation  (U)  is  satisfied  everywhere. 

The  system  of  differential  equations  (lU)  along  with  the  initial  conditions 
(12)  and  (20)  determines  a  coordinate  transformation 

(23)  X  -  x(t,cr)  ,  y  -  y(t,cJ)  . 

The  transformation   (23)  is  a  one-to-one  mapping  of  any  region   P  in  the   (t>  C) 
plane  into  a  region    P'      in  the    (x,y)  plane  provided  the  Jacobean  determinant 

is  different  from  zero  and  exists  throughout    P.     The   (tjCr)  coordinate  system 
is  orthogonal  as  one  can  see  easily  from   (21),    (lU)  and   (10).     The  rays  C  =  constant 
and  the  wavefronts    t  =  constant  of  the  geometrical  optics  system  fom  the  coor- 
dinate lines.     The  coordinate  system  is  singular  where  the  determinant  D  of 
(2li)  vanishes,   and  this  happens  generally  on  a  curve  which  is  the  envelope  of 
rays,   that  is,   a  caustic  of  the  ray  system. 

We  now  consider  the  geometrical  optics  amplitude  A  of  the  field.     Since  by 
(U)  the  magnitude  of  V^  is  n  we  can  write   (7)   as 

(25)  2n2  |A,   (^^)A  =0  , 

for  the  arclength  differential  ds  along  a  ray  is  given  by 

1 

ds  -  (X^2  .y^^)^di,~ld^, 
_ __ 

It  may  also  happen  on  a  focus  which  is  a  single  point  where  the  rays  intersect. 


according  to  (Ih),   and  the  dot  product  of  a  unit  vector  in  the  ray  direction 


dA 
3i- 


with  VA  is 

2    2  \ 

If  we  define  (Vtf)  =  v  and  take  into  account  (U)  we  can  express  the 

laplacian  A  in  the  ii,(i^)   coordinates; 

.    2  3^  ^  ,  2  d^  ^  1  d(n^D)  a  ^  1  a(v^)  d 

where  D  is  given  by  (2U).'-  -J     Thus  we  have 

(26)  At    .  I    ^J^    . 

After  a  substitution  from  (26)  into  (25)  there  results 


(27)  H-l 


»(^/(„2d)J*-0 


The  solirtion  of   (27)  is 

(28)  A  -  A^(cr)/(nD^^)  , 

where  A  (tf)  is  a  function  of  C  alone  and  is  determined  by  an  initial  condition 
for  A(i|r,<J')  on  each  ray  or  -  constant.  It  is  clear  from  (28)  that  A  beccwnes 
singular  on  a  caustic  where  D  vanishes. 


3.  The  Field  Amplitude. 

It  is  now  possible  to  add  to  the  differential  equations  (lU)  equations 
whose  solution  will  provide  the  quantity  D  of  (2U)  as  a  function  of  t  and  C,  so 
that  by  making  use  of  (28)  one  can  find  the  field  amplitude  A  along  any  ray» 


10  - 


If  the  rays  x(i|r,cr),  yi'ir,cr)   were  known  analytically  it  would  be  possible  to 
determine  D  directly  from  its  definition  (2U).  If  numerical  integration  methods 
are  to  be  used,  however,  as  matters  now  stand  D  can  only  be  estimated  in  this 
case  by  means  of  a  graphicail  representation  of  the  rays  and  an  appropriate  geo- 
metrical interpretation  of  D:  for  example,  in  its  role  as  the  ratio  of  crossec- 

[21 
tions  of  two  different  wavefronts  cut  out  by  an  infinitesimal  tube  of  rays.'--' 

Consider  a  quantity  M(t,cO  defined  by 
(29)  M  -  n^D  -  py^  -  qx^  , 

where  (2U)  and  (lU)  have  been  used.  According  to  (28)  ACt,©")  is  inversely  pro- 
portional to  the  square  root  of  M,  with  a  proportionality  factor  which  is  inde- 
pendent of  \|f  . 

Vath  the  help  of  (iLi)  and  (29)  the  derivative  M.  can  be  ejqpressed  as 

By  the  use  of  (29)  and  the  orthogonality  of  the  vectors  (p,q)  and  (x^,y^) 
(cf.  (18)  and  the  discussion  following  it]  x  end  y^  can  be  expressed  in  terms 
of  p,q  and  M.  That  is,  the  solution  of  the  linear  equations 

p^o-  -  qx^  -  M 

(31) 

qy^*  PV"  ° 

is 

x^  -  -  qM/n^ 
(32) 


V  "  ^/^^ 


11 


After  the  substitution  of  (32)  into  (30)  with  the  help  of  (lU)  there  results 


V  A  ^  p '  A  ^  M«  *  ^  si  '^/"''-  ^.  <"/"': 


The  derivatives  in  the  last  two  terms  on  the  right  of  this  equation  can  be 
expanded  by  the  use  of  (32)  and  the  last  equation  of  (lU): 

whence  from  the  last  equation  of  (lU)  there  results 


3(n^)  /t.2 


(33)  M^=j^^/(2n'q)JM-  p^q  . 

After  differentiating  the  second  equation  in  (lU)  with  respect  to  c  and 
making  use  of  (32)  we  obtain 


(3^>    v-Ap^[^/ng-jE-^/"5 


Equations  (33)  and  (3'!)  can  be  appended  to  the  system  (lU),  thus  enlarging  the 
system  of  ray  equations  so  that  a  complete  solution  will  give  p  and  M  as  well 
as  x,y,p  and  q  ^^ong  a  ray.  Initial  conditions  at  t  =  t  for  M  and  p  are 


(S^) 


M^(cr)  =  p^((?)yj((r)  -  q^(tf)x^(tf) 
Pcro  ■  Po(^)  • 
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Py  the  el imn nation  of  p  the  eq\iations  (33)  and  (3U)  can  be  reduced  to 
a  single  second  order  linear  differential  equation  for  H.  Thus,  after  solving 
(33)  for  p  and  differentiating  with  respect  to  ^   we  obtain 


yp^[^/n!l-^[^/4--n- 


(36)        %♦   3, 


From  the  substitution  of  (36)  into  (3I4)  there  results 


(37)  M   -  ^  A (log  n2)/(2n^)fM  -  0  . 

The  first  of  the  relations  in  (35)  provides  one  initial  condition  for  Vi   in  (37). 
The  other  initial  condition  can  be  obtained  from  (33)  and  the  second  equation 
of  (35).  The  initial  conditions  are,  then: 


^o^^^  '   Po^^^^o^""^  ~  %^^\^^^ 


(38) 


^(^)-ra/(-„\)k-p<; 


0/% 


where,    as  usual  the  subscript  o  attached  to  a  quantity  means  the  value  of  the 
quantity  at  t  =-  l'^. 

If  n     is  a  constant,    indicating  that  the  medium  is  homogeneous,    (37) 
reduces  to 


(39)  "^  "  °  • 

2 
The  same   equation  results  if  log  n     is  a  harmonic  function.     This  fact  is  con- 

2 
sistent  with  the  well-known  fact  that  when  log  n     is  a  harmonic  function  there 

exists  a  coordinate  transformation  under  which  the  wave  equation    (l)  for  an 


13 


inhomogeneous  medium  reduces  to  a  wave  equation  for  a  homogeneous  mediun,    that 
is,    vn.th  a   constant   index  of  refraction.      The  proof  is  simple j   for  the  sake 
of  completeness  it  will   be  given  here. 

Consider  the  coordinate  transformation 


(UO)  C  =  ^(x,y),       n=n(x,y) 

The  wave  equation   (1)   can  he  \irritten 


If 


(Ul)  (V^)^  =   (Vri)2   =  n^^      y^.VT]  =  C     and     A^  =  AT]  =  0  , 

the  V7ave  equation  ir;   (4,ti)   coordinates  becomes 

(U2)  ^^^  .^^^.k^,  =  0  . 

Now  the  conditions  (Ul)  can  be  satisfied  if  the  function  w(z)  =  4  +  in,  for 

z  =  X  +  iy,  is  an  analytic  function  of  z,  i.e.,  if  ^  and  tj  are  conjugate  harmonic 

functions  of  x  and  y,  and  if 

(U3)  Iw'(z)|2  =  n^  . 

2 

The  condition  (h3),  on  the  other  hand,  implies  that  log  n  is  the  real  part  of 

f  ,   1  2  2 

the  analytic  function  iv  {z}r     and*hence  that  log  n  is  harmonic.  The  steps 

in  the  proof  given  thus  far  are  reversible^  thus  it  follows  that  a  sufficient 

condition  for  the  existence  of  a  coordinate  transformation  (UO),  under  which  the 

wave  equation  (1)  for  an  inhomogeneous  medium  goes  into  the  wave  equation  (U2) 


lU  - 


2 

for  a  homogeneous  mediw.,  is  that  log  n  be  a  harmonic  function. 

To  explain  the  consistency  of  the  property  (39)  of  K  in  this  case  let 
I'l^  be  the  quantity  K  in  (29)  which  gives  the  field  amplitude,  according  to 
(28),  for  the  wave  equation  (Li2),  in  the  (4,ti)  coordinate  svstem,  and  let  M 
be  the  corresponding  quantity  for  the  wave  equation  (1)  in  the  (x,y)  coordin«te 
system.  Then  by  (2U)  and  (29)  we  have 


■     a(^,o-) 


and 

(U5)  M  -  n2  ^l3^  . 

a(^-,(3-) 

Now  the  well-lcnown  property  of  Jacobean  determinants  in  two  successive 
transformations : 

a(4,T)  ,  5(x,y)  ^  B(g,^) 
3(x,y)     d(!,<j)     d(n<y) 

implies  along  with   (bJj),    {h$),   the  Cauchy-Riemann  equations 

^x  "  V*       ^'-\' 
which  must  be  satisfied  when  4  and  T  are  conjugate  harmonic  funct-ions,   and  the 
conditions   (V^)^  «   (Vn)^  -  n^  that 

(U6)  M  -  K^   . 

According  tc  (U6)  the  quantity  K  must,  indeed,  satisfy  the  differential 
equation  (39)  if  >'  does,  end  M^  satisfies  (39)  because  the  medium  to  which  it 


-  13'  - 


corresponds  is  homogeneous,  with  index  of  refraction  one  according  to  (U2).  The 
above  discussion  proves,  incidentally,  that  the  condition 

2 
A log  n  =  0 

is  also  a  necessary  condition  that  there  exist  a  coordinate  transformation  of 
the  type  introduced  which  transforms  the  wave  equation  for  an  inhomogeneous 
medium  into  one  for  a  homogeneous  medium. 

Ii.  Examples. 

A,  The  Homogeneous  Kedium 

Suppose  the  medium  is  homogeneous  .so  that  (39)  holds.  Then  M  is  a  linear 
function  of  t  .  Thus,  by  (38)  we  have 


(U7)  M(1',(J)  -  -  lp^(o)/q^(o')^(t-t^)  +  M^(cr)  , 


where 


(U8)         M^(<r)  =  p^(cy)yj(cr)  -  q^(tf)xJ(o)  . 

According  to  the  standaini  geometrical  optics  result  the  field  amplitude 
on  a  wavefront  is  proportional  to  the  square  root  of  the  curvature  of  the  wave- 
front,  ^-^-^  By  comparison  with  (28)  and  (29)  this  has  the  implication  that  M  is 
proportional  to  the  radius  of  curvature  p  of  the  wavefront.  It  is  possible  to 
verify  that  (U?)  agrees  with  the  standaird  result,  then,  by  expressing  p  as  a 
function  of  ^  and  d. 

Let  t  be  the  arclength  parameter  along  a  ray  C  =  constant.  Then 

/,  o\       M    dM  /  dt   dM  /  ,  ,,  dx   ,,,  dy^   dM  /     dD 
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where  the  eiconal  equation  (h)  and  the  definition  (29)  of  M  have  been  used 
along  with  the  fact  that  V'/   has  the  direction  of  a  ray  and  the  fact  that  n  is 
constant. 

The  unit  vector  v  along  a  ray  normal  to  a  wavefront,  is  given  by 

(50)  V  -  (p,q)/n  . 

If  we  differentiate  v  with  respect  to  the  arclength  parameter  s  along  the  wave- 
front,  we  obtain 

«^'      £  -  i  S  l§  (1.  -  P/<"  -  K  H  {(H)'  *  (g' j  'ft.  -  P/^'  • 

A  result  of  differential  geometry  is:     the  magnitude  of  -r-  is  the  curvature,   or 
the  reciprocal  of   the  radius  of  curvature,    of  the  wavefi-ont.     Thus,    after  making 
use  of    (32)  and  the  last  equation  of   (8),   we  obtain  for  the  radius   of  curvature  p 
of  the  wavefront 

(52)  p=    l(Mq)/(n|£)|    =    |(nl>l)/|£   |    , 

where  the  definition  (29)  has  also  been  applied.  Since  p  is  the  distance  from 
the  point  of  tangency  of  a  ray  on  a  caustic  to  the  point  of  intersection  of  the 
wavefront  with  the  ray,  p  differs  fr^m  the  arclength  t  by  a  constant,  and  thus 
dp  -  dt.  On  the  other  hand,  since  n  is  constant,  from  (lU)  it  follows  that 
p  and  q  are  independent  of  ^  ,    i.e.,  functions  of  0"  alone.  These  facts  and 
(U9)  and  (52)  imply  that 

(53)  M^-tpJ/q^. 


17 


The  sign  in  (53)  can  be  deterirjined  by  observing  that  -r—  has  the  direction  of 
the  vector  (x  >y_J>  but  according  to  (32)  the  vector  M(q,  -  p)  has  the  opposite 
direction.  From  (51)  and  (52)  it  follows  then  that  (53)  must  be 

which  verifies  the  result  (U?). 


B •  The  Medium  Whose  Index  of  Refraction  is  a  Function  of  One  Variable. 

2 

Suppose  that  n  is  a  function  of  x  alone  and  that  a  plane  wave  is  incident 

2  2 

from  -  00 J  i.e.,  n  will  be  required  to  approach  a  constant,  n  and  the  solution 

|j,  of  (1)  approaches  exp  j  ikn  (x  cos  a  +  y  sin  a)!-  as  x  approaches  -  oo .  In 

2 
addition,  it  will  be  assumed  that  n  is  smooth  and  everywhere  positive  and  that 

2 

lim   ''j  ^   =  0  .  The  initial  curve  will  be  taken  to  be 


(5U)  y  =  cr,  X  =  -  00. 

The  quantities  p  and  q  will  be  constant  initially: 

(55)  p(-  00 )  =  n  cos  a,  q(-  oo  )  »  n  sin  a  . 

Fran   (II4)  it  follows  that  q  is  independent  of   if  and  thus  will  have  the 
same  value  evei'j'^here.     It  then  follows  from  the  last  equation  of   (II4)  that 


(56)  p  -  \/n^  -  nf  sin^a 


Equation  (37)  can  be  expressed  in  terms  of  the  variables  x  and  y.  If  it 
is  assumed  that  M  is  independent  of  y,  the  resulting  equation  will  be  consider- 


Ifl  - 


ably  simpler  than  it  would  be  without  this  assumption.  The  final  result  will 
justify  this  assumption  because  of  the  uniqueness  of  the  solution  of  (37),  (38). 
From  the  first  equation  in  (lU)  and  equation  (56)  there  results 

(57)  M^  .  x^M^  -  pM^n^  -  |n^  -  n^^  sin2a|  Mx/n^ 
After  using  (57)  in  (37)  we  obtain 

(58)  (n^  -  n^  sin2G)(MVn^)'  *  }^\r?)' /{2r?)   -  (n2)7n2[  M  -  0 

where  the  primes  stand  for  differentiation  with  respect  to  x.  Equation  (58) 
can  be  rewritten  as 


(59)  1(n^  -  n^sin2a)M'/n^f  -  \  {n^ )\/ (.2r? ))     -  0  . 


After  one  integration  (59)  becomes 

(60)  (n^  -  n^^  s±r?a.)n  /n     -   (n^)  V(2n^)  -  B  , 
where  B  is  a  constant.  The  general  solution  of  (60)  is 

1   X  ,3  1 

(61)  M  -  B(n^  -  n^^sin^a)^  J  n^  (n^  -  n^^  ain^a)  ^dx  +  C(n^  -  n^  sin^a)^  , 


sith  B  and  C  arbitrary  constants.     From  the  initial  conditions    (38),   the  initial 
surface  givBn  by   (5U)  and  equations   (55)  there  results 


(62)  M(-  00  )  -  n     cos  a  ,     M  (-  oo  )   =»  -  3-  cot  a  ■  0  . 

o  '        \|f  dy 
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The  first  equation  of  (62)  and  equation  (61)  imply  that 

(63)  C  -  1  . 

2 

The  second  equation  of  (62)  and  equation  (61)  imply  since  lim    '  "  ^  =  0, 

X  ->-00 

that 

(6U)  B  =  0  . 

Thus  we  have,  finally, 

1 
(65)  M  =  (n^  -  n^  sin^a)^  . 

Using  (28)  and  (29)  we  obtain  the  field  amplitude: 

1  _  1 

(660  A  -  (n^  cos  o)^(n^  -  n^  sin^a)  ^  , 

where  the  first  equation  of  (62)  and  the  fact  that  the  field  amplitude  approaches 

one  as  x  approaches  -  oo  have  also  been  used.  The  result  (66)  agrees  with  that 

[o] 
of  Seckler  and  Keller •- J  who  calculated  the  geometrical  optics  solution  of  this 

problem  using  the  standard  method  for  obtaining  the  field  amplitude,  that  is, 

computing  the  square  root  of  the  ratio  of  wavefront  cross  sections  cut  out  by 

an  infinitesimal  tube  of  rays  on  the  initial  and  final  wavefronts. 
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